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EXPLICIT CONSTRUCTION OF NON-STATIONARY FRAMES FOR L 2 

UBERTINO BATTISTI, MICHELE BERRA 


Abstract. We show the existence of a family of frames of L 2 (R) which depend on a parameter 
ae [0,1]. If a = 0, we recover the usual Gabor frame, if a = 1 we obtain a frame system which 
is closely related to the so called DOST basis, first introduced in m and then analyzed in [4]. If 
a s (0,1), the frame system is associated to a so called ^-partitioning of the frequency domain. 
Restricting to the case a = 1, we provide a truly ^.-dimensional version of the DOST basis and 
an associated frame of L 2 (M d ). 


1. Introduction 


One of the most intriguing problems of modern Time-Frequency analysis is the construction of 
new efficient methods to represent signals, which can be one dimensional or, more often, multidi¬ 
mensional, such as digital images. The increasing amount of data and their complexity forces to 
develop optimized techniques that address the representation in a fast and efficient way. 

The starting point of this paper is the definition of the S-transform, first introduced by R. G. 
Stockwell et al. in [3Tj as 


( 1 . 1 ) 


(S/)(ft,0 = (2«r i Je- 2 " it£ |e|/(t)e- 


r(t-br 


dt. 


The expression of the S-transform (1.1) is formally very similar to the so called Short Time Fourier 
Transform or Gabor Transform (STFT) with Gaussian window 


(STFT /) (b, 0 = (27t)-* J e- 2 ^f{t)e~^dt. 


The main novelty of the S-transform is the frequency depending window. The leading idea is the 
heuristic fact that, in order to detect high frequencies, it is enough to consider a shorter time. 
Therefore, the width of the Gaussian is not fixed but depends on the frequency, shirking as far as 
the frequency increases. The S-transform was introduced to improve the analysis of seismic imaging, 
and it is now considered an important tool in geophysics, see [I]. In |2Ij . the connection between 
the phase of S-transform and the instantaneous frequency - useful in several applications - has been 
studied. See 0 © m m on luma na for some applications of the S-transform to signal processing 
in general. 

From the mathematical point of view, M. W. Wong and H. Zhu in [35] introduced a generalized 
version of the S-transform as follows 

(S* /) (6,0 = f e- 2 ™* / (t) |f | <p{£{t-b)) dt, b, £ e R, 

Jr 


where ip is a general window function in L 2 (R). The S-transform has a strong similarity with the 
STFT, actually it is also possible to show a deep link with the wavelet transform, see El 1321- 
In fact, the S-transform can be seen as an hybrid between the STFT and the wavelet transform. 
Representation Theory provides a very deep connection among S-transform, STFT and wavelet 
transform, as they all relate to the representation of the so-called affine Weyl-Heisenberg group, 
studied in [23]. This connection has been highlighted in the multi-dimensional case by L. Riba in 
[25) . see also [25]. The affine Weyl-Heisenberg group is also the key to represent the a-modulation 
groups, [7]. 
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Our analysis focuses on the DOST (Discrete Orthonormal Stockwell Transform), a discretization 
of the S-Transform, first introduced by R. G. Stockwell in 30]. In [3], the DOST transform has 
been studied from a mathematical point of view. It is shown that the DOST is essentially the 
decomposition of a periodic signal / G L 2 ([0,1]) in a suitable orthonormal basis defined as 

(1-2) |J D P = U > D r>-r}r=0 > 

peZ peZ 


where 


2 P — 1 

Dp, T {t) = J—- Yi e 27UTI ( t ' \ t G R,p ^ l,r = 0,... ,2 P_1 - 1, 

v2 p v= 2 p-i 

with the convention that Do(t) = 1, see Section [3] for the precise description of the basis func¬ 
tions. The DOST basis has a non stationary time-frequency localization, roughly speaking the time 
localization increases as the frequency increases, while the frequency localization decreases as the 
frequency increases. Therefore, the basis decomposition of a periodic signal is able to localize high 
frequencies, for example spikes. The time-frequency localization properties of the DOST basis imply 
that the coefficients 

/p,r = (/> ^P,t)l 2 (\o,i\) > 

represent the time-frequency content of the signal / in a certain time-frequency box, which is related 
to a dyadic decomposition in the frequency domain. Moreover, this decomposition can be seen as 
a sampling of a generalized S-transform with a particular analyzing window which is essentially a 
box car window in the frequency domain. 

The DOST transform gained interest in the applied world after the FFT-fast algorithm discovered 
by Y. Wang ang J. Orchard, see [33j 134] for the original algorithm and [4] for a slightly different 
approach which shows that the fast algorithm is essentially a clever application of Plancharel The¬ 
orem. 

The orthogonality property is clearly very useful for several applications, nevertheless it is well 
known that, in order to describe signals, a certain amount of redundancy can be very effective. 
Therefore, it is a natural task to look for frames associated to the DOST basis. 

We follow the idea of the construction of Gabor frames: 

G(g,a,/3) = {T ak M 0n g} = I g (t - a k)\ 

^ J ( k,n ) eZ 2 


If g is a suitable window function, for example a Gaussian, and a/3 < 1 then Q(g,a,/3) is a frame, 
see for example [2D]- It is possible to consider the Gabor frame as the frame associated to the 
standard Fourier basis, which is formed by all modulations with integer frequency, extended by 
periodicity and then localized using the translation of the analyzing window function g. Inspired 
by this approach, we consider the system of functions 


(1.3) [T k ^(D p , 0 (t)g(t))}^ k ^ 22 , P(p) = 2^ -1 if p ^ 0, /3(0) = 1, i/ > 0. 


For simplicity, here we have not introduced a frequency parameter. The idea of the system in (1.3) is 


to consider the DOST basis D p o and then localize it with a window function g. The main difference 
from the standard Gabor system is that the translation parameter is not uniform, but depends 
on the frequency parameter p. Therefore, (1.3) can be considered as a non stationary Gabor frame, 


using the terminology introduced in Ej. Inspired by the theory of a-modulation frames, see El. 
El. El. El; in Section [ 3 ] we introduce a family of bases of L 2 ([0,1]) depending on a parameter 
a G [0,1], If a = 0 we recover the standard Fourier basis, if a = 1 the DOST basis; when a G (0,1) 
we show that the basis is associated to a suitable a-partitioning of the frequency domain in the 
sense of [ID] . 

In Section |4j we prove the main result, see Theorem 4.13 we show that for each a G [0,1] the 
localization procedure explained above produces frames of L 2 (IR), provided the time and frequency 
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parameter are small enough. The main tool is a non stationary version of the Walnut representation, 
see Subsection 14.21 

In Section [5j we analyze the higher dimensional case. Restricting to a = 1, we consider a 
multidimensional partitioning of the frequency domain and the associated frames. This construction 
is different to the usual extension of the DOST to higher dimensions. In |33j|34|, the DOST applied 
to two dimensional signals (digital images) was essentially the one dimensional DOST applied in 
the vertical and then in the horizontal direction, therefore it was not a truly bi-dimensional version 
of the transform. 


2. Notations 


In the paper we use the following normalization of the Fourier Transform 

F(/)(ca) = /M = J e~ 2ni ™f(x)dx. 

In the multidimensional case, we do not write explicitly the inner product in R d , we leave the same 
notation of the one dimensional case. 

Set 

M k f(x) = e 2nikx f(x), k,xeR d . 

and 

T k f(x) = f(x-k), k,xe R d . 

Let f,geL 2 (R d ), we denote the L 2 -scalar product as 

</) 9> = J fgdx. 

We denote S(R d ) the Schwartz space. 


(3.1) 


3. a-BASES of L 2 ([0,1]) 

Let a e [0,1], we define a partition of N associated to the parameter a. We use an iterative 
approach. Let p be a non negative integer. For each a, we define 

^a; 0 0, Sa;0 1; /^a(0) 1, 

^a;p S<x;p— 1? ^a;p ^a;p T [^ a; p] Pa ( P ) lda;p-l’ P ^ 

where [xj is the integer part of x. Then, set 

I a - p = [*a;p, Sa;p) n N pth interval of the partition of N, 

Pa (p) = \Ia;p\ = |so;;p — *a;p| width of the interval I p . 

We then set for all p e N 

-i Saw— I 

„ m =_ 

P,T 


B*{t) = 


XI e 2mi, ( f t = 0 ,..., p a (p) 


1. 


VA*(p) 

See Figure [l] for a plot of real and imaginary part of such basis with different values of a and p. 
For p e Z_ we define 


= B_p, T> T = 0, . . .,Pa(-p) - 1. 

Notice that if a = 0 then /3°(p) = 1, for all p , therefore r is always zero, hence 

B%(t) = e 2 ™ pt . 

That is (J Bp 0 is the ordinary Fourier basis of L 2 ([0,1]). 

If a = 1, then ii- 0 = 0, Si ; o = 1 and then i\. p = 2 P_1 , si ;p = 2 P , Pi(p) = 2 P_1 forp ^ 1. Therefore 

-®o,oW = 1 
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Figure 1. Elements of the DOST bases. Notice that a = 0 is the classical Fourier 
basis which has no localization in time. As a or p increases we gain time localiza¬ 
tion. The black line and the dotted-red one represents the real and imaginary part 
respectively. 


while 


2 P —1 


<r(t) - t 


T = 0, 


2 p - 1 - 1 . 


r/=2 p ~ 1 


That is (J p T Bp T is the so called DOST basis, see [5], [501 . 

For now on, we restrict to positive integers, all the results hold true also for negative integers via 
simple arguments. 

Clearly the above partitioning can be considered on R instead of Z. Let us consider the parti¬ 
tioning of the real line. Notice that for a = 0, trivially 


0°(p) = 1 = V°, Vije/o-.p- 
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Figure 2. Growth of f} a ( p ) with the respect of the parameter a. 


If a £ (0,1] then f} a (p) = |i“ p J and for r\ e I a . p = [i a;p , i a - p + |*“ p J) we have, for |p| large enough, 

1 „ L*2;pJ ^ l*2;pj „ *2;p _ 

2“ + 1 " (2 i a ., p T " ( V) a " *g;p ' 

Using Fornasier’s notation, see [16], we can write 

|la;p| ~ \v\ ! V ^ Icr,p- 

That is the above partitioning is an a-covering. 

Remark 3.1. Notice that the a-partitioning introduced above is well defined for all a £ [0,1]. The 
case a = 1 is not defined as a limit case, as in the usual analysis of a-modulation spaces, see mm- 
The key point is that we use an iterative scheme, instead of a definition involving the function . 
In this way, we can get rid of the singularity which arises at a = 1. At a = 1 the growth of fi(p) 
with respect to p is exponential while for a £ [0,1) is polynomial, see Figure [£| Using the iterative 
definition of (3{p) this fact causes no problems. 

Theorem 3.2. The functions 


R“ t (-), p £ Z, t = 0,..., /3 a (p) — 1 
form an orthonormal basis of L 2 ([0,1]). 

Proof. Notice that for a = 0 and a = 1 the Theorem holds true in view of well known properties of 
Fourier basis and of results proven in 3]. 

The proof follows closely the argument in 0]. 

Step 1 I ^p,r | l 2 ([ 0 , !]) = 














6 


UBERTINO BATTISTI, MICHELE BERRA 


Since {e 2mkt } fe( _ z is an orthonormal basis of L 2 ([0,1]) we can write 




i 2 ([o,i]) p t 


i s “iP 1 Sa 'P 1 r 

33 J J J 


3 2 ™ raw) e" 2 ™ ^'C*- fcTpi) dt 


1 


ri'=i a - p 

Sa;c 1 Sq; tj 1 




1 


^q: ;z> 1 


A* (p) 


S 1 


= 1. 


Step 2 (Bp T , B C *, T ,) 


L 2 ([0,1]) 


= 6 0 (p-p') So (t - r'). 


If p = p' and t = t' the first step implies the assertion. Moreover, if p ^ p' then I a . p and I a . p i are 
disjoint, thus 


(i?“t,£ p V) 


L 2 ([0,1]) 


= 0. 


So, we can restrict to the case p = p'. Since {e 27Tlfc *} feeZ is an orthonormal basis of l? ([0,1]), we 
obtain 


( B p,T, B p,T') L2m]) ^ ^ 2 




Pa(p) 


rj=i 

Pa (p)-l 


/3a (?) 


^ e 27ti(r—r') 


(*Q!;p+j) 
0a-(?) 


i=o 


Let us suppose r — r' # 0, then we can write 

e 2 ™( T - T ')Krifr 

/ TDOi TDOi \ _ __ 

V p// L 2 ([0,1]) 


Pa (?) 


^ ^ e 27ri(r-r') 


/\ l \ 3 


Pol (p) 


j=o 


(3.2) 


s 2Tri(r-T')^\ / x _ e 27Ti (t-t') 

Jll- e 2Tzi(T ~ T ' ] ^ 


Pa (?) 


= 0. 


In equation (3.21 we used well known properties of geometric series and the fact that J ^ is not 


an integer number and therefore e 2 ™( T-T )/Pa(p) ^ \ 
Step 3 The functions B p T generate L 2 ([0,1]). 
Notice that 


Pa (p)-i 

U - s P an { e2 ™ fct }fee[i Q;p ,s Q; p-l] 


r=0 


Therefore, to prove that the functions U“ T are a basis of L 2 ([0,1]) it is sufficient to check that 
3 a 

Pv 


Uriif) 1 33p T arc linear independent. That is, if 
(3.3) 


Pa (p) 1 

2 OL T B p T = 0 => a T = 0, r = 0,..., pa (p) - 1. 


T = 0 
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Let us consider the projection of (3.3) into the Fourier basis. We obtain the system of equation 
P a (p )-1 


(3.4) 


z 

T=0 


OL'r & 


— 27ii r 


Mp) =0, j = 0,.. .,p a (p) - 1. 


We can rewrite the above equation as a linear system 

e - 2 7ti(/3 ct (p)- 1 )^2J 

— 27Ti - 


(3.5) 


(1 
1 


— 27li a a } P ' 
e Pot{P) 

P+1 

Pot (p) 




V i 


l '.\pP~Pot (P) — 1 
Pot (p) 


= -27Ti(/3 Q (p)-l)- 


t\p+Pot (p) — 1 
Pot(p) 



( «o \ 


( 0 \ 


Ot\ 


0 


V a /3c(p)-l ) 


l 0 ) 


The square matrix in (3.5) is a Vandermonde matrix. Since the entries (e 2m ( l ^,p+j)/P Q ( p ))^.“ 0 


&<x{p)~ i 


are all distinct the determinant of the matrix is zero and therefore the system in (3.51 has only the 
null solution, that is in (3.4) and (3.3), a T mus vanish for all r = 0,..., /3„ ( p ) — 1. □ 


3.1. Localization properties of the functions H“ T . Now, we investigate the time-frequency 
localization properties of the functions B“ T . They clearly have compact support in the frequency 
domain and the support is precisely I a - p . Therefore, they cannot have compact support in the time 
domain. Nevertheless, it is possible to determine localization properties in the spirit of Donoho-Stark 
Theorem, P33- 

Property 3 . 3 . For each a e [ 0 , 1 ] and p , t, the following inequality holds 


T _i_ 1 

’Po(p) + 2 / 3 ct (p) 


I B^ T \ Z dt 


^ 0.85. 


\ Pot(p) 2Pot (p) 

that is the L 2 -norm is concentrated in the interval 

T 1 


J3 a (p) 2/3 a (p) ’ 13 a (p) 2f3 a (p) 

if r = 0, the interval must be considered as an interval in the circle: 


0, 


1 


2/3a ( P ) 


u 


1 - 


1 


2/3a ( P ) : 


Proof. The proof is based on Taylor expansion and Gauss summation formula. In [J the Property 
is proven for the case a = 1, actually the same proof works also for general ae [0,1]. □ 

4. DOST- WAVE PACKET FRAME, DIMENSION d = 1. 

In order to give a more flexible structure to the DOST basis, we generalize it to a redundant and 
non-orthogonal system and show that this leads to a L 2 -frame. First, we address the 1-dimensional 
case in full generality taking into account the phase-space tiling dependent from the parameter a 
defined in the first part of the paper. 

4.1. Frame definition. The idea is to extend by periodicity the DOST basis and then localize 
them using a suitable window function. 

Definition 4.1. Consider 


(4.1) 


= 


1 


nr -Z e 2 ™^ (p,k) e Z 2 , 

v p a {p)j 


where I a - P is defined in Section [3| Then, we define the a-DOST system 

D a (p, v, tp) = {<Pp, k (t)} Pik , (P, k) e Z 2 , 

see Figure^} or the plot of the frame elements. We did not plot the case a = 0 which is the standard 
Gabor frame element. 
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Remark 4.2. Notice that the Fourier transform of ( |4.1| ) simplifies into 

1 —27ti ujk 




V P <*(p) 


“ / 3 c* (p) z (p{oj - rj) 

p€!ll a . p 


( P ,k)eZ 2 . 


This kind of system of functions has affinities with the well known family of Non-Stationary Gabor 
frame, see HIED H2]. 

The function 

M= X! V{u-rj) 

r]€pl a - p 

mimic a bell function of the set I a . p as shown in Figure [^J 

Remark 4 . 3 . Sometimes, when we want to highlight the dependence of the element ip^ k (x) on the 
frame parameters, we use the notation Tf,’k’ V • 

Remark 4 . 4 . As we pointed out is Section [3J when a = 0 the DOST basis reduces to the Fourier 
one. As one may expect, when we analyze the same case for the a-DOST system above, we end up 
with a Gabor frame. Precisely: 

D ° (/bG T>) = G(p, G P) = {TpkMppip, ( p , k) e Z 2 }. 

Indeed, recalling that (3 a ( p ) = 1 when a = 0, i a . p = p we thus obtain 

<*(*) = Z e 2 ^ x - vk \{x - uk) = T vk M pp p. 

V=PP 


We investigate the condition under which the family D a (p., v, p) is a frame of L 2 


4.2. Walnut-like Formula. We retrieve a Walnut-like formula for our DOST system in the spirit 
of Gabor analysis. This representation is a very useful tool to prove the frame property. Moreover, 
when a = 0, we recover the Walnut representation for Gabor frames. 

First, we recall the Poisson formula, J2U] Proposition 1.4.2]. 

Lemma 4.5. Suppose that for some e > 0 and C > 0 we have \p(x)\ C(1 + |x|) _d_£ and 

|<£(w)| C(l+ \u\)- d ~ e . Then 

(4.2) 2 <p(x + 7") = 2 ^ (-) ^ nX, \ 

neZ d neZj d ' ^' 

for 7 > 0 , where d is the dimension of the base space. 

We can state the main result of this section. 


Lemma 4.6. Let tp,if e L 2 (M.), f e Cf? ( K), and 
(4.3) 


p,k 


s^f) (t) = z</> t e k, (p, k ) e z 2 . 

(?) 


Then 

(4-4) F t _» w (w) = ^ Z Z (/ M) 

where 


peZ fce 1 


i’p-.a [ U) - k- 




VEUla-,p 
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(a) a = 0.5, p = 5 



(c) o = l, p = 5 



(b) a = 0.5, p = 7 



(d) q=1, p = 7 


Figure 3. Plot of the elements of the a-DOST frame in time, with a Gaussian 
function as a window and [i = 1/2. As we did for the bases, we observe that the 
localization increases when a or p grows. The black line and the dotted-red one 
represents the real and imaginary part respectively. 


Proof. Take the Fourier transform of the frame operator and obtain 


m= Z a<r>cr m 


p,ke Z 2 


V /f(.) 1 M\ 1 

^ VPAp) 


e 27 tife 3^o“^. a (w) 


p,ke Z 2 

=zfz</(-). 


1 k , ^ . M ,, / w 1 27ii k a s uj 


'J Pa ip) 


e 27tife ti.np)M^ ;a (.)) 


VMp) 




= " _1 Z 7TTT Z</(0.^ fc ^ ( - ) ^-(0>e anifc ^“ 

V / 
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(a) a = 0.5, p = 5 (b) a = 0.5, p = 7 




(c) a = 1 p = 5 (d) o=l, p = 7 

Figure 4. Plot of the elements of the DOST-frame in frequency, with a Gaussian 
function as a window and p = 1/2, v = 1/5. Notice that as p grows, the window 
translates and approximates a bell function of the set J a;p . 


Then, (4.2) with p = /<f> p;Q ,, n = k and 7 = yields 

F x ^ {S$$ v f) M = (w - k^-) tt£. a (w) 

peZIceZ V ' ' 

= S S / (" -*—) *?»("- ) ®s»m 

peZ IceZ ' ' ' ' 

= z'- 1 E E (/>)) ^ (w - $j a (w), 

tip?, fce7. ^ \ / 


pe Z /cgZ 


as desired. 


□ 


Remark 4.7. ITe stress that the equalities (4.3) and (4.4) are, at this stage, formal. The frame 
property we prove in Theorem f.ll implies that, under suitable assumption on the function p, 
is a linear operator from L 2 (R) to itself, therefore (4.3) and (4.4) are equalities in L 2 (R). 


4.3. Painless frame expansion. In order to have a clear scheme of the frame construction we 
first consider the painless case. 


Definition 4.8. Let p e <S(R), such that supp</ c \—L, L] and 0 < p < 1. Given p e R, p > 0, we 
say that (p, p) is admissible if 
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( 4 -5) ^;«H= Yi 

satisfies the following properties 

1. 0 < %. a {uS) Ci, for some Ci > 0. 

2. For every uel \{p : 4>^ ;q ,(u;) > 0}| < C 2 . 

3. There exists C 3 > 0 such that for every uj e R d , $^ ;a (w) > C 3 , for some p e Z. 
Using this definition we can immediately obtain some important properties of . 

Lemma 4.9. Let (<p,p) he admissible, then 

(4.6) 

(4.7) 


|supp <f>£. Q I (2L + p) p a (p ), Vp e Z, 
C3<I]l^;aMr<C?C'2 ) V W G 


where C\,C- 2 ,C 3 are defined above (cf. Definition ^, tty . 

Proof. The first property comes directly from the definition of (cf. ( |4.5| )) and the fact that 
supp<p cz \—L,L\. First, we observe that ui — p G suppp if — L ^ u> — r) ^ L; thus, writing 
rj = 1 u, ( i a - p + t), this translates as 

pia;p — L ^ UI < /rs a; p + L. 

Hence 

|supp<F£. Q | < p (s a - p - i a . p ) + 2 L = p/3 a {p) + 2 L. 

Since L is fixed and /3 a (p) ^ 1, uniformly in p, we can write 

| supp $p ;a | <(p + 2 L)/3 a (p ), 

as desired. 

Given ui G R, we know that there are at most C 2 indexes for which 4>((. a (w) > 0. Moreover, 
^ Ci, thus the upper bound of (4.7) holds true. 

On the other hand, for the same u>, there exists at least one p such that $p ;a (w) ^ C 3 . □ 

Remark 4.10. An explicit example of an admissible function is the compact version of a Gaussian; 
precisely <p(ui) = Xe( a ; )^2 e_7IW > w ^ ere Xe a bounded smooth function such that 

/ \ f 1) W 6 [ — 1; 1] ! 

XeM = n r 1 ■, , 1 

^ 0 , id E [— 1 — £, 1 + • 

In |27j . similar functions are considered. 

Theorem 4.11. Consider (p, p) being admissible. Then, there exist positive lattice parameters v,p 
such that the a-DOST system (cf. (4.1)^ 

D a ( p , v ,< p ) = {< fe (*)} (pjfc)eZ2 

is a frame for L 2 (R). Precisely, there exist A, B > 0 such that for all f G L 2 (R) 

(4.8) m\ 2 2<Y\<f^pA x ))\ 2 ^ B \\ni 

p,k 


Proof. The frame property (4.8) is equivalent to 
Hence, we evaluate 




<Wf,f> = <FWf),f>, 
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where F is the Fourier transform as above. The Walnut representation formula (4.4) yields 
<F (5“;^/) ,/> = ^ T 

peZ feeZ 


(4.9) 


^ pez 

+ ^<£ £ T ^P)f T iMp)^K; a Jy 


peZfceZ\{0} 


We notice that by (4.7) 


Indeed, 




;<£ I $ p;«| 2 />/> " </./> = II/II2- 


We claim that for v small enough, 

1 


-<£ £ T $M P )f T ZM P )^K,<*’f> = 0 - 


peZ fceZ\{0} 


Recall that |supp4>((. a (w)| < (2 L + fj,)/3 a ( p ) (cf. (4.6)), thus for v < (2L + [T) 1 and k ^ 0 


suppT^^^a n supp$£ ;a = 0, 


as desired. 


□ 


4.4. Conjugate Filter. We can define a way to represent functions even without a canonical dual 
frame; namely, we construct a conjugate filter for the window function ip. This technique is a 
powerful tool for numerical implementations, see 123- 
Set 


(4.10) 


= * 


$p;aM 


£ feM 


12 ’ 


then, the product between and forms almost a partition of unity 
(4.11) = 


We have the following corollary: 


Corollary 4 . 12 . Consider the functions defined in (14.10h , then set 


'J'p.fcM = 


1 


—27ti cuk 


V/3a(p) 


Qp£(oj) 


Then, for any f e L 2 
(4.12) 


/W = £</.^, fc K,fe(t)- 

p,k 
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Proof. We notice immediately that supp = supp (w). We take the Fourier transform of 

(4.12), thus 

E</> *£*>%(*) 

p,k 


2 (2 (a 


(P) > </,e “"‘TOT] (felpj-te )*»(«) 


peZ \fceZ 

(4.13) = 2 (^(p)- 1 2 F (/«£«) (fc. 


\ — 7T —2niujk 7 . 


keZ 


fy “ ™ K>)- 


Then, by Poisson formula (cf. (4.2)), (4.13) turns into 


2 (""‘I;(/« 

peZ \ fceZ 


w + fc 


/3a (p) 




Repeating the same procedure as in the proof of Theorem 4.11 i.e. excluding the terms with k > 0, 
we can conclude that, for v small enough, 


(4.14) 


E</,*“ fc >0“ fc (*) = -/M E n ^M^ ; a(w)- 

p,/c peZ 


By (4.11), equation (4.14) becomes 


p,k 


Finally, applying the inverse Fourier transform, we obtain (4.12). 


□ 


4.5. DOST frames, general construction. In this section we prove that we can build up a a- 
DOST frame with milder hypothesis on the window function compared with the ones of the previous 
section. 


Theorem 4.13. Consider a function ip e L 2 ( R) n L 1 ( R). As above, set $^. a (w) as (cf. ( |4.5| ) ) 

= E 

pe/i/ Q; p 

Suppose that 

a<E l^;aMI 2 <&, V.eK, 

P 

for suitable constants a,b> 0. Moreover, assume that the window ip satisfies 

(4-15) IftMl < (i + C |^|)iv ’ VweR 

for some N > 2. 

Then, there exists A, B > 0 and > 0, such that for all f e L 2 (R) 


m\u\<s^f,f>\^B\\f\\ 


2 

25 


for 0 < v < vo • 
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4.6. Preparatory Lemmata. We need some result concerning the decay of the elements 4 + a 
outside the sets I a - P . 

Lemma 4.14. Let ip e L 2 (R) n L 1 (R) such that: 


I^HI < 


C 


TV 


(i + M) N ’ 


with N > 2. Then, 
(4.16) 


\KM\ < 


c 


N,/i 


(1 + d (v,fil a . p )) 


TV—1 - 


Proof. From the definition of we obtain: 


C 


N 


Hence, 


E E 7, ,, _ ,,AT 

rje/ a ;p J7e/ Q ; P "b | w t* 7 ?!) 


I^Mlsd + dK^r” 2 c K (yd( u .ui » 

r/el a;p "b l w TV\) 


If w G pI a - P , then d (w, pI a -, P ) = 0. Thus 
(4.17) 


<*>£»!< £ cw(i + ]w-Hr*<avE(/M + i) 

vel a;p je Z 


-A 


which is clearly convergent for N > 2. 

If ui f pl a ;p, then d (w, pl a ;p) = 6 and 


I J| ^ (1 + 6)N £ Q (l + rt + 6)N 


< 


2Cn 
(:1 + S) N 

2Cn 

~ (1 + 6 + 

2 Cn 

~ (1 + 6 + 

2Cn 

~ (1 + 5 + 

for IV > 2, which completes the proof. 
Lemma 4.15. Let wet and N > 1. Then 


1 + 

1 + 

1 + 


(1 + 6 ) 


TV 


0 (1 + (IX + b) N 

+ oo ^ 


dx 


o (1+T&*)" 

+ 00 


dx 


(1 + 6 ) 


f 

JO 


T Jo (! + x) N 

Cn,p 


dx 


i + M C ]< 


(1 + 6 ) 


TV-1 


E 


nTV 


^ C 


pez (1 + d (u,fxl a -p)) 
and the constant is independent on u. 

Proof. For any wel, there exists only one p such that ui G pl a -,p- Then, for j ^ 2, 


3~ 1 


d (w, tll a] p+j) > J]/3 a {p+ k) ^ (j 


□ 


fc=l 


1 ), 
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since (3 a ( p ) ^ 1. The same argument works for j < —2. 
Obviously 


V - 1 - 

7=— 1 (1 + d (w, pl a -p + j)) 


^ 3. 


Hence, we can bound our sum as 


£ 


i 

(1 + d (w,pI a;p )) N 


y Cn_ 

Zj aN 

jeZ\{0} J 


and the constants do not depend on the particular choice of w, as desired. 


□ 


Lemma 4.16. Let as in Theorem 4-13 above, then 


^ £ E 


4>p;a ( OJ - k 


Pa(p) 




peZ feeZ\{0} 

Proof. Let us define for all p e Z 


= 0. 


and 


We can write 


(4.18) 

(4.19) 

(4.20) 

(4.21) 




E £ 

peZ IceZ\{0} 


4>p; a IW-fc 


&*(p) 


<£ £ 

peZ feeZ\{0} 


ldlpja I UJ k 


Pa{p) 


+ £ £ 

peZ fceZ\{0} 

+ £ £ 

peZ fceZ\{0} 

+ £ £ 

peZ fceZ\{0} 


2 < l ) p;a ( /t 


I 4>p;c 


2 4)^ ;q , ( w - k- 


*(p) 


* 0 ) 


Ap) 


l*£;a M 

2^„M 

2^;aM 


00 


For now on, let us suppose vp. < 1; then we notice that the term in (4.18) is identically zero for 
each p e Z, since the supports are disjoint. Clearly this is not restrictive since we are considering 
the limit for v —* 0. We remark that: 

(4.22) d - k ^ a ^ P \ p,I a -pj > - pj /3 a ( p ), Vw 6 pl a - p , 


with 


|fc| — up > 0, fee Z\{0}. 











































16 


UBERTINO BATTISTI, MICHELE BERRA 


We then analyze the term in (4.19). For each ui there exists a unique p = p such that ui G pl a; p, 
hence for each ui 


z z 

peZ feeZ\{0} 


2^»p;a W-fc 


.fa) 


iKa M 


- z 

fceZ\{0} 


$?.„ [id — k 


Pa{p) 


l$5;u M 


^ C 


N,n 


z 

fceZ\{0} 


C N 


,p 


( 1+d ( 


(4.23) 


< 


< 


Cn,„ 2 


C N 


N—l 


feeZ\{0} (l + | (i£I - /i) Pc 


N—l 


C N ^C N ((Z _ ^ 

fceZ\{0} ' ^ 


1-AT 


We made use of the decay proven in Lemma 4.14 and of (4.22). Since N — 1 > 1 and the inequality 

>1 


(4.23) does not depend on w, we obtain 

Pa(p) 


z z 

peZ feeZ\{0} 


2^p;a 


l^p;a M 


(4.24) 


< Cn,ijC n 2 

k€ Z\{0} ' ^ 


= C n ^C n v n 1 2 


1—N 


1 


fceZ\{0} (1^1 


1—N ' 


Since the sum is convergent uniformly with respect to u, the term in ( 4.24| ) goes to 0 as v —* 0 with 
the rate v N ~ x . 

In order to consider the term in (4.20), notice that for all w G R, there exists a unique p such 
that id G pl a -p. Moreover, arguing as above, for all p G Z\ {p} there exist a unique k p G Z\ {0} such 
that to — k p g I a - p . Hence, we can write (4.20) as 


(4.25) 


z 

peZ\{p} 


By (4.22) 


1 4>p; Q I Id — k 


\k 


Paip) 


2$£;aM 


d(id, pilot- p) ^- p) Pa (p), if Id -k. 


Pa ip) 


G pla\p- 


Finally, by (|4.16|), we can write 
«(P)' 


' / co P el\{p} [ 


peZ\{p} 


r 12 


(l + ^w.^Jaa ,))^- 1 




z 


U N,p, 


peZ\{p} 




(4.26) 


z 

z 


(1 + d(id,fiI a .p)) N 1 e (1 + d(td,pl a -p)) e 

C n,p 1 


peZ\{p} + ^(^••/' / n;p)) V -' * (l + (M - ^ & 

C 


peZ\{p} 


(1 + d(id,pl a -p)) 


N-l-e ' 


The sum in (4.26) converges uniformly with respect to id by Lemma 4.15 therefore the whole sum 
goes to zero as v does. 
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The last term to be taken into account is (|4.21|), which includes the “tails” of our window function. 

7^ e ///«;?, the 

2^;a(w) = 0. 


Fix peZ, by definition if either to e pI a - tP or w- e pl a . p , then 

0 a {p) 


Define the set 


2' 1 ' P ;a 


= U : 0 < d 


^Z/cT a;p j < ^/3 q (p) . 


Notice that only a finite number of element k belong to this set; precisely, when pv < 1, then 

uniformly in ui and p. 


K\ < 2, 


So, for all to, we can split the term in (4.21) as follows: 

Pa(p) 


z z 

peZ fceZ\{0} 


2 '*p,a 




^Z z 

peZ k$F“ 

+z z 

peZ fceFy 


2 < l > p;a I OJ A: 


*(p) 


2 ^p;a 


0a{p) 


2$£aM 




We show that both terms are bounded by a constant times v e with e > small enough independent 
on 10 . 


We notice that if k £ F p , then 


, /3a (p) T 

IO - k - ,pl a ;p 


\k — k\ . . 

> J -/3a (p) 


where k is the closest index to k inside F£\ hence, we can rearrange the (infinite) indexes k e Z\{0} 
into j e Z\{0} such that 


■ 0a (p) 


d W - j - -,/x/a;p ^ — 0 a (p) ■ 


\j\ 


From the discussion above, and the estimate (|4.16|) it follows that 

0 a(p) 


z z 

peZ fc£Fy 


2 %-a W - k 




^Z z 

peZ jel\{0} 


c 


N,p 


c 


N,fi 


^1 + ^0a (p)'j (1 + d (lO, pl a -,p)) 

Cn,p Cn,p 


< 


< 


N—l 


y y 

peZ jeZ\{0} (M) " 1 (1 + d (w, pla-p)) 


N—l 


z 


Z 


je Z\{0} 


N — l 


C N ,p 


3 / (1 + d(w,p/ Q;p )) 


N -1 


The latter term is summable in j, Lemma 4.15 implies the summability in p as well, therefore the 
whole term goes to zero as v N ~ l . 

The last part follows from the observation below: if there exists k e F£, then 


0a ip) 


1 


0 < d - k - — — , pla\pJ < 7 ^ 0 a ip) ■ 

Since p is fixed and v —» 0, we can assume pv < 1; thus, there exists 5 such that pv = 1 — 6. We 
split our analysis in two cases, first we consider 
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Let s e fj,l a;p such that d (u>, s) = d (w, /i/ a;p ) then by triangular inequality 


d (u>, pl a ;p) = d (w, s) ^ d I W, LO — k 


Pa ( P) 


Pa ( P ) 


>\*\ 


Pa ( P ) 


— + P) Pa(p) . 


Hence 


A , T Pa(p) (i £>\ paip) (\ n 1 , 5 

d (u},pl a . p ) > \k\— -— ( 1 — 2 ) = M fc l - 1 + 2 


If |fc| O 1, it is clear that 


d (u, fil a . p ) > —p a (p) 


Then, 




c 


N,fi 


Cn, p.v' 


(1 + d (u,pI a;p )) N 1 e ^i+| PsMY (1 + d (u,pl a .p)) 
On the other hand, when 

, ( ,Pa{p) T \ § o , x 

d I W - k —-— ,pl a . p 1 > — Pa (p) 


N-l-e ■ 


we have 


2^p;ck 


«(p) 


< 


Cn^ 


< 


(l + d (w- k^P,pI a -^ 
Cn,p 


N -1 


which goes to zero as . Using the inequalities, it is clear that in both cases 


(4.27) 


E E 

peZ fceF w 


2 ^p;o: 


*(p) 


,3+ 

^ p;a 


M 


< 


+“i 2 


C N v e 


pez (i + d (w,/x/ a;p )) 


N—l—e 


since the sum is uniformly bounded with respect to w, the above inequalities imply that (4.27) goes 
to zero as u e . The quantity (1 + d (w, /r/ a;p )) ~ 1 ~ £ is summable if N — 1 — e > 1. This is granted 
by the fact that N > 2 and that we can chose e small enough. 


□ 


4.7. Proof of the main result. 
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Proof. As we did in the painless case, we can represent the frame operator as (cf. (4.91) 

tip) 


/> = < F (wf), ;> = -<e E T k m (/ m) 

pe Z fceZ 

= M,/M> 


( 00 — k~ 






+ ^<E E 

^ peZ fceZ\{0} “ \ v / 

HEI*&« Ml 2 />)./»> 


^p;o; 


< 


+ ;2 2 

peZfc e z\{ 0 } 


X 2 


* 0 ) 




<T fc ^) (/m) ,/H) 


peZ fceZ\{0} 


4>p;a I 00 — k- 


M 




T !; Mp) 


-7. 2 


peZ IceZ\{0} 


4>p;q ( w - k- 


M 


(/ m) 

*£;«M 


m. 


Similarly, 


<s$pf,r» ‘l, i 5f2l^;«(“)r-2 2 

peZ feZ\{0} 


1/ ueR 
V p 

By hypothesis, we know that 




M 




ll/ll 


Lemma [4.16| implies that 


lim V' V 

iy->0 ^ ^ 


pe Z 

Pa{p) 


K<* [ U - k~ 


peZ fceZ\{0} 

Hence, there exists vq > 0 such that for all 0 < v < io 0 

M 


E E 

peZ feeZ\{0} 


4>p;a ( W — fc- 


$p;aM 


KM 


= 0. 


<: 


Then, the action of the frame operator can be bounded as follows 


l<s“r/>/>!<HEIfeMllL + E E 

IpeZ peZ fceZ\{0} 

< l ( 6 + I) H/ll 2 


<T>P 

p,ot 


M 


KM 


ll/ll; 


l<5;,r/./>|»Hi^2l <|, t«MI 2 -2 2 

peZ fceZ\{0} 


$p;a I W - k 


PM 


KM) 


> & 11 / 112 - 


as desired. 


□ 
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4.8. Existence of DOST frames. We show that the Gaussian satisfies the hypotheses of the main 
Theorem 14.131 


Theorem 4.17. 


Consider p = e 


and set $p ;a (w) 


as in (4.5). 


Then, 


(4.28) 


weR, 

V 


for suitable constants a,b > 0. Moreover, 

(4.29) I^MI< (i+^tivv , 

for any N ^ 0. 


Proof. The polynomial decay claimed in (4.29) is trivial, since ip is a Schwartz function. For the 
lower bound in (4.28), we argue as follows: for any wet there exists only one p such that ui e pI a ; P - 
For the sake of simplicity, we assume p > 0, the negative case follows with the same argument. 
Therefore, 

P 

Since uj e then uj = + £), 0 ^ t ^ /3 a (p). Hence 

2 


l$- (« 

I T)\Ot \ 


W)l =2 

_ 1 
“ 2 

_ 1 
“ 2 


E 


,-n^n(i a . p +t)-ri) 7 


r/efj,I a; p 


da (p) — 1 

^ e -nWap,p+t)-n(.i a -,f+j)y 

3=0 

da (p)-l 

E e 

3=0 


2 


1 

> - 
2 

sup e ^ 

3 = o, —iA»(p)-l 


because of the positivity of the Gaussian. The maximum value is reached when 1 1 — j\ is small. Our 
construction implies that there exists j such that \t — j\ ^ 1, thus 


which is independent on p and u, as desired. 

Due to the fact that the Gaussian is positive, it follows that also is positive as well. Hence 



P P 

We rewrite the sum above as follows: 

p j 

Finally, since ip belongs to the Wiener Space, we can write 

( \ 2 2 

esssup^]<£(a;-/xj)J < ((~ + IMIwJ < °°> 

where || • ||^ is the Wiener norm. We have used a well known property of Wiener space, see e.g. 
|20| [Lemma 6.1.2], 

□ 
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5. Higher Dimensions 

We consider here the case a = 1 and an arbitrary dimension. We define a (parabolic) phase 
space tiling and, for suitable window functions, we provide a frame of L 2 (M d ). We follow the ideas 
of wave atoms proposed in M and subsequently adapted to the Gaussian case by E23- For the 
sake of simplicity, we enlighten the notation used before by suppressing the parameter a. 

As for the dimension d = 1 , we begin with the painless case using a smooth and compactly 
supported window function. Moreover, we can define an explicit conjugate frame that leads to a 
reconstruction formula. Then we generalize the construction. 


5.1. Phase space partition. Define the Cartesian coronae C p as follows: 

Co = [—1, l] d , 

C p = \uj = (wi,... ,u d ) G : max |w s | e [/3 (p), /3 (p + 1))) , p ^ 1. 


l^i^d 


Each corona is further partitioned in (open) boxes of side /? (p) = 2 V 1 , precisely 

d 

X p/ = Y[[f3(p)-£ S: p(p)-(e s + l)), 


s =1 


where i s = — 2 , — 1 , 0,1 and max s= i..i ( l s + ^) f Co, i.e. the centers are outside the inner corona. 
The indexes l s label every possible box inside the corona. It can be easily checked that the number 
of such boxes (or multi-indexes) is 2 d (2 d — 1), for every p ^ 1. We also define, according to Section 

0 

(5.1) 


Ip = x P ,t f> d , 


Figure [5] 


see 

We generalize now the 1-dimensional DOST system (cf. (4.1)). 

Definition 5.1. Given the set I^ (cf. & , consider 

(5-2) ipyix) = 0(p) - 2 V e 3 ™ , K*~*5fio) v > (x- kfC] , 

fin v 

where 7 = (p, k, £) e (N, Z d , J), J contains the admissible indexes £ described above. Then, we define 
the (multi-dimensional) DOST-system 


MD (p, v, tp) = {<p T (x)} 7 . 

See Figure [fi| as an example of frame element both in time and frequency. Figure [?] shows how the 
elements are localized in the set I d in frequency. 


5.2. Painless frame expansion. As for the case of dimension d = 1, we start with compactly 
supported window functions. We adapt the definition of admissible window. 

Definition 5.2. Let p e <S(R d ) such that supp<p cz [— L,L\ d , for some p > 0 and 0 ^ p < 1. We 
say that (p, p) is admissible if, given 


(5.3) 


= E - V) 


it satisfies the following properties 

1. 0 < 4>p.^(u;) ^ Ci, for some Ci > 0 . 

2. For every u<=R d \{p,£: $^(w) > 0}| < C 2 . 

3. There exists C 3 > 0 such that for every uj e R d , 4C.^w) > C 3 , for some (p, £) e N x J. 
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P ip + 1) 


Pip) - 


i e 

P 


- + 


X 


p,i 


p ip) 


Pip+1 ) 


Figure 5. 2D phase space partitioning, and the set 


Using this definition we can immediately obtain some important properties of 
Lemma 5.3. Let i<p,p) be admissible, then 

(5.4) Diam (supp sc Cl^P ip) , Vp e N, 


(5.5) 


p,£ 


where C\, C%, C 3 are defined above (cf. f.8) and Diam denotes the maximum distance between points 
of the set, i.e. the diameter. 

Lemma 5.4. Let ip,if in L 2 (M d ). For each f e L 2 (K d ), set 

(S%%f) W = 7 e (N,Z d ,J). 


Then 

(5.6) 

where 


S' 




p,£ k 


K* {u ~ k 


Paip) 




K,-A u} ) = Z 


The proofs of these lemmata are analogous to Lemmata |4.6||4.9| 

Theorem 5.5. Consider (<p,p) being admissible. Then there exists v > 0 such that the DOST- 
system 

MD (p,u, ip) 

is a frame for L 2 (M d ). Precisely, there exist A, B > 0 such that for all f e L 2 (R d ) 

(5.7) m\i^J]\<f,pAx)>\ 2 ^m\i 

7 

The proof is analogous to the one made in dimension d = 1. 
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5.3. Conjugate Filter. Set 

(5.8) 

then 


= " d 


KM 


S \KM 


(5.9) 


Y>KrMKM = » d - 

p,l 

Corollary 5.6. Consider the functions Kf'f defined above, set 


Vh 7 ( W ) = 0(p)~* e- 2mu *»=T!l*r(uh 7 = (p,k,£) e (N,Z d , J), w e : 


Then, for any f e L 
(5.10) 


7v--y r- w - --p; 

2 /in>d\ 


/(*) = 


Proof. We briefly review the proof of Corollary |4.12| Nothing changes here a part from the normal¬ 
ization constant. 

X!</’ 

7 

= E (WH </H,e- 27tiwfc ^>^M>) ^(p)“5 e- 27ti “ fc ^)$^H^ 

= 4E ( E e - 27riw<; TO^ 


\keZ d 


Then, by Poisson formula (cf. (4.21) and ( |5.9| ), supposing u 0 small enough, we have for all 0 < v < u 0 

p,£ 

as desired. □ 

5.4. DOST frames, general construction. We state and prove that we can build up a DOST 
frame with similar hypothesis to the ones given in Section [4] 

Theorem 5.7. Consider a function ip e L 2 (M. d ) n L 1 (K d ). Let 3>p.^(w) (cf. (5.3),) 

KM = E M~V)- 

pep-Tp 

Suppose that p > 0 is chosen so that 

a^J]\KM \ 2 VcueR, 


for suitable constants a, b > 0. Moreover, assume that the window tp satisfies 

( 5 .ii) l«“)l < 

for some N > 2d. 

Then there exists A, B > 0 and isq > 0, such that for all f e L 2 (K d ) 


m\ 2 2<\<s^fj)\^B\\f\\\ 


for some 0 < v < v$. 
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5.5. Preparatory Lemmata. We recall the same result proved in dimension d = 1. 


Lemma 5.8. Let ip e L 2 (R d ) n L 1 (R d ) such that (5.11) holds true. Then 

C N 


(5.12) 


I^MI < 


(l + d(uj,nl*)) 


N-d ' 


Lemma 5.9. Let w e and N > d. Then 


£ 


i 


P ,e (1 + d (w,/*/*)) 

and the constant is independent on w. 


N 


^ c 


Lemma 5.10. Let '3?^ as before then 


S 


P,i feeZ d \{0} 




= 0 . 


Proof. As we did in Lemma 4.16 set 


MM = XplMKM 


and 


2 KM = <!>£.» - 

Then, we split the norm as follows: 


I £ 

P,i fceZ d \{0} 




KM 


(5.13) 

(5.14) 

(5.15) 

(5.16) 


= E £ 

P,* fceZ d \{0} 


Md ( w - fc - 


: (p) 


+£ £ 

fceZ d \{0} 

+ £ £ 

P^ fceZ d \{0} 

+ £ £ 

p,e feeZ d \{0} 


2 ^p;< ( W-fc- 


1 $^ ( u - k- 


M* I w 


Hp) 

V 

] (p) 

V 

>(p) 


MM 


MM 


2 KM 


2 KM 


Notice that if vp, < 1, then the term in (5.13) is identically zero for each p £ N, since the supports 
are disjoint. 

In order to analyze the term in (5.14), notice that for each uj £ R d there exist unique p,l such 
that oj e iL Notice that we have 


(5.17) 


k/3 ( p ) 


> \k\ Q - Vdfi) /3(p), Vw 


£ I. 

p 


For each ui £ M , by (5.17) we have 


M-.M 2K 


p\£ 


' - M) 


< 


r 2 


(l + d (w, pip)) 


N-d 


i — VdfJ^j |fc|^ 


d—N 
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Our hypotheses grant that N — d > d, then the above remarks implies that 


E E 

P,l /ceZ d \{0} 


i ®p-,e M 2 $^ luj-k 


Hp) 


< 


1 2 ^[u-k 


'(p) 


E x/di - 1 ) k 


fee Z d \{0} 


d—N 


< 


C Q - Vd/r) 


d-M 


where the latter tends to 0 as u —> 0 and the constants are uniformly bounded with the respect of 
to. Hence the term in (5.14) has a limit vanishing as v approaches zero. 

The term in (5.15), goes to zero as v goes to zero as well. We notice that for each u>, there exist 
p, l such that ui e i|. For the same reason for each p, l e (N\ {p}, J\ {F}), there exists k p j such that 
w — fcFM 6 H Therefore, 


(5.18) 
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Then, applying equation (5.17), (5.12) and Lemma 5.9 as in (4.26) we can conclude that the term 
in (5.18) goes to zero as v goes to zero independently on u>. 

We consider now the term in (5.16), set 


Fp f = -(fc: 0 <d|od — 


kfi (p) 


,P T p) < i^P(p) 


then 
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uniformly in u>. 


One can split the term in (|5.16|) as follows: 
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Then, for each ui 6 . 


and 


This yields 
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which goes to zero as v goes to zero as desired. We stress that (l + d (u,pIp)) N d 6 is summable 
if N — d — e > d which is granted by Lemma 5.9 Since the bounds are all unifrom with the respect 
to w, we can conclude that the terms in (5.16) goes to zero as v does. □ 
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5.6. Proof of the main result. 

Proof. From the d = 1 case, we get 

< l { sup H 0 (u) + H/III 

(s^f, /> > h {jg, Ho(w) - H" k> 1 {uj) } ll/ll! 

with 

H o = J]\K-A “)\ 2 
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By hypothesis, we know that 
Lemma 


5.10 


a ^ Hq{lS) ^ b. 

implies that lim 1/ _ >0 #£>1 = 0- Hence, there exists v 0 > 0 such that 


i < f, VO < i/ < Z2 0 . 

Then, for all ^ e (0, i/q), the action of the frame operator can be bounded as follows 


\<S^ v f,f)\ < i {||HoL + II/II2 < l(b+ |) It/ll 




//n - II. 


1 


II/II2 > ^ ll/ll 2 


2 ' 


□ 


Remark 5.11. As we did in dimension d = 1 (cf. 
Gaussian fulfills the hypothesis of Theorem\5. ?\ 


Theorem f.ll), we can show that the normalized 


6. Conclusions 

We constructed a frequency-adapted frame which covers Gabor and Stockwell-related frames. 
Our setting includes also general o-phase-space partitioning. This approach appears natural to 
describe a-Modulation spaces and this will be subject of a future work. 

In [3j , the author prove that the DOST basis is able to diagonalize the S'-transform with a suit¬ 
able window function which is essentially a boxcar window in the frequency domain and that the 
evaluation of the DOST-coefhcients turns out to be the evaluation of the S'-transform with this 
particular window in a suitable lattice. The natural question is if the a-DOST bases introduced in 
this paper have the same property, clearly not with respect to the S-transform, but in relation to 
another transform, similar to the flexible Gabor-wavelet transform (or a-transform ), see e.g. (16l 
for the the definition. 

The n-dimensional case considered in Section [5] is restricted to the case a = 1, hence a suitable 
phase-space partitioning is yet to be defined for a e [0,1) and will be part of our future studies. 
This issue has been already analyzed in the two dimensional case by N. Morten in [25] using the 
theory of Brushlets. 

From a computational stand point, we aim to implement and compare our results with existent 
methods. We are interested in testing in various applications such as medical and seismic imaging 
and also general image processing. As pointed out in the introduction, we remark that our approach 
consider the 71 -dimensional case in a peculiar way: instead of applying the one dimensional DOST 
in each direction sequentially, we provide a native 71 -dimensional setting. This approach is similar 
to the Wavepackets and Curvelets one, see mm 

A natural question arises: is the density of our frames comparable with the Gabor case? For 
instance, is it true that if the volume of the lattice is strictly lower than 1, the Gaussian leads to a 
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Figure 6 . Time and frequency outlook of two window functions. We observe 
heavy decay in time while in frequency we localize around a certain frequency. The 
window tp is in both cases a normalized Gaussian. 


frame? And is this condition independent on a? 
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